August,  1976 


Report  ESL-R-683 

Contract  ONR/N00014-75-C-1183 


SHORTEST  ROUTE  ALGORITHMS  FOR  SPARSELY 
CONNECTED  NETWORKS 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY,  Cambridge,  Massachusetts  02u» 

Department  of  Electrical  Engineering  and  Computer  Science 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dal.  Enle,.*;  UNUljAbS.lt1  IED  


•*  * ► REPORT  DOCUMENTATION  PAGE  befoIeDcomple™form 


1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


TITLE  (and  Subtittm)  

'SHORTEST  .ROUTE  ALGORITHMS  FOR  SPARSELY 
CONNECTED^  NETWORKS  , ~ ' 


7.  AUTHORfal 


5.  TYPE. OF  B6PQRT  ft  PERIOD  COVERED 

Technical j 


6.  PERFORMING  ORG.  REPORT  NUMBER 


5) Joe  E.  A)efenderfer 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Massachusetts  Institute  of  Technology 
Electronic  Systems  Laboratory 
Cambridge,  Massachusetts  02139 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Defense  Advanced  Research  Projects  Agency  1 

1400  Wilson  Boulevard 

Arlington,  Virginia  22209  ' 


1«.  MONITORING  AGENCY  NAME  ft  ADDRESS  (II  dllloront  from  Controlling  Olllco) 

Office  of  tlaval  Research 
Information  Systems  Program 
Code  437 

Arlington,  Virginia  22217 


16.  DISTRIBUTION  STATEMENT  (ol  thl,  Report) 

Approved  for  public  release;  distribution  unlimited. 


K UNIT  NUMBERS 

Program  Code  No.  5T10 
ONR  Identifying  No.  049-383 


\Sm.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


0&ESL-T-4>f3 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


19.  KEY  WORDS  ( Continue  on  reverse  elde  If  necessery  end  identify  by  block  number) 


Shortest  routes 
Networks 

Decomposition  algorithms 


20  ABSTRACT  (Continue  on  rave raa  tide  II  necesaafy  ,nd  Identity  by  block  numbmr) 

This  report  studies  the  shortest  route  problem  for  networks  that  are  less 
fully  connected.  Two  algorithms  are  presented  which  exploit  the  absence  of  arcs 
in  solving  the  shortest  route  problem.  The  first,  which  is  desiqnated  the  NXN 
algorithm,  would  tend  to  be  the  more  applicable  to  networks  typically  encountered 
in  practice.  The  second,  which  is  an  improvement  on  Hu's  decomposition  shortest 
route  algorithm,  is  more  efficient  for  a small  class  of  networks;  however,  it 
generally  requires  less  memory  to  hold  the  required  decomposition  information  in 
the  computer  than  does  the  NXN  algorithm.//  


August,  1976 


Report  ESL—R— 683 


SHORTEST  ROUTE  ALGORITHMS  FOR 
SPARSELY  CONNECTED  NETWORKS 
by 

Joe  E.  Defender fer 


Q □ 


This  report  studies  the  shortest  route  problem  for  networks  that 
are  less  than  fully  connected.  Two  algorithms  are  presented  which 
exploit  the  absence  of  arcs  in  solving  the  shortest  route  problem. 

The  first,  which  is  designated  the  NXN  algorithm,  would  tend  to  be 
the  more  applicable  to  networks  typically  encountered  in  practice. 

The  second,  which  is  an  improvement  on  Hu's  decomposition  shortest 
route  algorithm,  is  more  efficient  for  a small  class  of  networks; 
however,  it  generally  requires  less  memory  to  hold  the  required 
decomposition  information  in  the  computer  than  does  the  NXN  algorithm. 


Acknowledgements 


I wish  to  thank  Professors  John  M.  Wozencraft  and  Robert  G.  Gallager 


who  provided  encouragement  and  many  helpful  discussions 


This  work  was  supported  by  a Sloan  Research  Traineeship  and  the 


Advanced  Research  Projects  Agency  under  Contract  ONR/N00014-75-C-1183 


r 

! 


..7. 


I 


-4- 


Table  of  Contents 


Abstract 

Acknowledgements 


Page 

2 


Section  I 

Introduction 

7 

II 

The  NXN  Algorithm 

9 

III 

Decomposing  the  Network  for  the 
Algorithm 

NXN 

12 

IV 

The  IHU  Algorithm 

15 

. 

V 

Decomposing  the  Network  for  the 
Algorithm 

IHU 

20 

1 

VI 

Some  Examples  Using  the  IHU  and 
Algorithms 

NXN 

21 

« 

Appendix 

' 

• 1 

• 

Section  IA 

Introduction 

27 

| 

IIA 

Bookkeeping  for  Shortest  Routes 

28 

■ 

IIIA 

The  Main  Program 

29  ' 

IVA 

Subroutine  DIJKST 

32 

f 

VA 

Subroutines  DECIHU  and  IHU 

33 

i 

VIA 

Subroutines  DECNXN  and  NXN 

37 

» 

VILA 

Program  Listing 

42 

Biblography 

71 

1 

Distribution  List 


72 


5- 


Figure  1 
2 

3 

4 

5 

6 

7 

8 

9 


An  N node  network  with  an  NXN  decomposition. 

The  network  of  Figure  1 with  a distinct  NXN 
decomposition. 

The  form  of  the  D matrix  for  the  IHU  algorithm. 

The  topology  of  ARPA  network  with  decomposition 
information. 

A 47  node  symmetric  network  with  decomposition 
information. 

A 64  node  network  with  decomposition  information. 

An  example  network  with  seven  nodes  and  thirteen  arcs. 

Node  renumbering  and  partitioning  by  subroutine  DECNXN 
for  the  network  of  Figure  7. 

The  topology  of  Figure  7 with  new  node  numbers  as  assigned 
by  DECNXN. 


Tables 


Table  1 


Comparative  performance  of  three  different  shortest 
route  algorithms  on  the  three  sample  networks. 


Topology  cards  for  the  network  of  Figure  7. 

Cards  punched  by  subroutine  DECIHU  for  the  network 
of  Figure  7. 

Values  of  variables  in  common  block  IHUSTF  deduced 
from  Table  3. 

Cards  punched  by  subroutine  DECNXN  for  the  network  of 
Figure  7. 

A Fortran  Program  which  interprets  the  cards  punched 
by  DECNXN. 


-7- 


Section  I Introduction 


The  problem  of  finding  all  the  shortest  routes  in  a directed  network 


has  an  extensive  literature  [3,9]  due  to  the  number  of  network 
problems  to  which  shortest  route  algorithms  are  applied.  This  paper 
presents  two  new  shortest  route  algorithms  which  can  significantly  reduce 


the  required  computation  time  when  the  network  is  less  than  fully  connected. 
The  first  is  based  on  original  decomposition  ideas  and  is  called  the  node- 
by-node  decomposition  (NXN)  algorithm.^  The  second  is  based  on  Hu's  de- 


composition algorithm  [5,6,11]  and  is  designated  the  improved  Hu  (IHU) 
algorithm. 

The  shortest  route  problem  is  formulated  as  a shortest  distance  problem 


where  D = [d^  ] is  a given  matrix.  The  number  d^  represents  the  length 

of  the  directed  arc  from  node  i to  node  j,  and  thus  it  is  assumed  d = 0. 

A path  P from  i to  j is  an  ordered  sequence  i=k,k, ,...,  k , ,k  = j, 

o 1 m-1  m 

m 

and  the  length  of  the  path,  L(P) , is  defined  as  L(P)  = £ d If  P 

r=l  r-1  r 

is  any  closed  path,  then  it  is  assumed  L(P)  > 0 so  that  the  shortest  distance 

problem  is  well  defined.  Then  the  problem  is  to  find  D*  = [d* . ] where 

^ ] 

d?_.  = min  L(P)  for  P ranging  over  all  paths  from  i to  j.  Knowinq  D*  alone 
does  not  specify  the  shortest  routes,  but  is  a well  documented  fact  that 
by  appropriate  bookkr  eping  as  one  calculates  D*,  the  shortest  routes  can  also 
be  established. 


After  completion  of  this  paper,  the  equivalence  of  the  NXN  algorithm  with 
previous  work  done  at  Network  Analysis  Corporation  under  ARPA  Order  No.  1523 
[8]  was  discovered. 
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Typically  D*  is  calculated  as  a series  of  refinements  on  D.  Floyd's 
algorithm  [4]  is  cited  for  an  N node  networks 

For  every  i e{l,2, . . . ,N>,  do  step  a: 

a)  For  every  j,k  e{1,2, . . . ,n} , do  step  bs 
b)  djk  - min  <djk,  dji  + d.k) 


where  " means  "is  replaced  by".  The  algorithm  requires  N additions 
and  N3  comparisons,  and  it  is  generally  assumed  additions  and  comparisons 
take  about  the  same  amount  of  time  so  that  one  says  Floyd's  algorithms 
requires  2N3  operations.  At  the  conclusion  of  the  algorithm  D*  has  replaced 
D.  Proof  of  the  algorithm  is  found  elsewhere  [6],  but  the  interested  reader 
can  easily  convince  himself  that  when  i has  been  stepped  from  1 though  i 

o 

then  the  current  value  of  d^k  is  the  minimal  distance  over  all  paths  from 

j to  k under  the  condition  that  the  intermediate  nodes  are  elements  of  the 

set  {1,2,...,  l }. 

o 

Ho  alqorithm  which  solves  the  shortest  route  alqorithm  could  be  any 

simpler  to  encode,  but  there  are  a variety  of  faster  algorithms  in  terms 

of  number  of  operations  [7,10].  The  standard  aqainst  which  the  new 

decomposition  algorithms  will  be  measured  is  Yen's  implementation  of 

3 3 

Dijkstra's  alqorithm  [2,11]  requiring  operations.  The  algorithms 
claiming  even  less  operations  are  not  significantly  faster,  theoretically, 
for  networks  of  the  size  for  which  computational  experience  is  cited  in 
this  paper?  furthermore,  some  of  the  apparent  gains  of  the  theoretically 
faster  algorithms  would  be  offset  by  their  additional  algorithmic  complexity. 


Section  II  The  NXN  Algorithm 


The  NXN  algorithm  for  solving  the  shortest  route  problem  is  actually 
a special  case  of  the  following  new  2N3  operation  algorithm: 


1)  For  every  i£{l,2, . . . , N-2}  in  order,  do  step  a: 

a)  For  every  j,ke{i+l,i+2,...,N},  do  step  b: 

b)  d.,  min  (d.,  ,d..+d  ) 

3k  3k  lk 

2)  For  every  ie{N-2,H-3, . . . ,l)  in  order,  do  step  a: 

a)  For  every  j ,ke(i+l, i+2, . . . ,n} , do  steps  b and  c: 

b)  d . . min  (d . , +d,  . , d . . ) 

13  lk  k]  13 

c)  d . . min  (d  . , +d,  . , d.  . ) 

31  3k  ki  31 

An  intuitive  proof  of  this  algorithm  will  be  helpful  in  understanding  the 
NXN  algorithm.  By  inductive  reasoning  similar  to  that  for  Floyd's  algo- 
rithm, when  step  1 has  been  completed  for  i=iQ,  then  d (for  j,k  > iQ) 
represents  the  conditional  shortest  j to  k distance  subject  to  all  inter- 
mediate nodes  beinq  elements  of  the  set  {l,2,...,i  }.  Consequently,  when 

o 

step  1 has  been  completed,  then  the  d (for  j,k  > N-2)  represent  uncon- 

ditional  shortest  distances  d* . . 

13 

Note  than  an  arbitrary  i to  j path  (for  j > i)  must  be  of  the  form 
i, . . . ,r , . . . , j where  r is  the  first  element  in  the  path  such  that  r > i; 
and,  if  this  path  is  the  shortest  path,  then  its  length  is  d|r+d* . When 
performing  step  2 for  i = N-2,  d*_.  is  known  and  d?^  must  be  the  same  as 


the  set  {l,2,..., N- 3 } ; it  follows  that  at  the  end  of  step  2 for  i = N-2, 

d,  . = d*  , and  similarly  d..  - d* . for  every  je{N-l,N}.  Clearly,  inductive 
13  i}  31  31 

reasoning  shows  that  at  the  end  of  the  algorithm  D = D*. 

The  NXN  algorithm  will  now  be  presented.  However,  in  order  to  simplify 

the  discussion,  it  is  assumed  that  all  of  the  arcs  are  duplex,  i.e.  if  d^<°° 

then  d..<°°.  Define  C.,  called  the  ith  connection  set,  as  follows:  jeC.  if 

31  l l 

j > i and  there  exists  a path  P from  i to  j such  that  L(P)  < 00  and  every 
intermediate  node  k satisfies  k < i.  Notice  that  the  are  functions  of 
topology  only  (implicitly  assuming  the  length  assigned  to  an  arc  is  00  if 
and  only  if  the  arc  does  not  exist  in  some  sense) . 

In  step  1 of  the  above  algorithm,  d..  = 00  if  j0C.  and  d = 00  if 
k^C^.  Furthermore,  in  step  2 of  the  above  algorithm,  = 00  and  d^  = 00 

if  k0C^.  The  corresponding  operations  are  clearly  unnecessary;  the  algo- 
rithm obtained  by  deleting  them  is  called  the  NXN  algorithm: 


1)  For  every  ie{l , 2, . . . ,N-2)  in  order,  do  step  a: 

a)  For  every  j,keC^,  do  step  b: 

b)  d *■  nin  (d . , d..+d  ) 

3k  3k  31  lk 

2)  For  every  ie{N-2,  N-3,...,l}  in  order,  do  step  a: 

a)  For  every  jc{i+l, i+2, . . . ,n}  and  k£C^ , do  steps  b and  c: 

b)  d.  . *■  min  (d  +d  . ,d..) 

ij  ik  k3  13 

c)  d . . *■  min  (d  ..  +d.  . ,d  . . ) 

31  3k  ki  31 
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A decomposition  is  defined  as  an  ordering  of  the  nodes.  Since  the 
connection  sets  are  a function  of  the  decomposition#  the  number  of 
operations  which  the  algorithm  requires  is  also  a function  of  the  decom- 
position, as  will  be  demostrated  in  the  following  section. 

In  the  case  where  some  of  the  arcs  are  not  duplex,  two  alternatives 


are  available.  The  first  is  to  change  the  definition  of  Ch  as  follows: 


jeCh  if  j > i and  there  exists  a path  P from  i to  j or  from  j to  i such 


that  L(P)  < 00  and  every  intermediate  node  k satisfies  k < i.  This  approach 
causes  unnecessary  operations  for  the  algorithm.  The  alternative  is  to 
define  two  connection  sets  for  each  node — one  for  the  incoming  connections 
and  one  for  the  outgoing  connections.  In  the  latter  case,  one  must  alter 
the  NXN  algorithm  to  incorporate  the  efficiencies  of  the  additional  connection 
sets.  The  increased  algorithmic  complexity  of  the  second  approach  and  the 
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Section  III  Decomposing  the  Network  for  the  NXN  Algorithm 

This  section  is  introduced  via  an  example.  Consider  figures  1 and  2 
in  which  the  same  network  has  been  decomposed  two  ways.  For  the  first, 

Ci  = {i+l,N-l,N}  when  ie{l,2, . . . ,N-3>  and  CN_2  = {n-1,n}j  the  number  of 
operations  for  the  NXN  algorithm  is  calculated  in  a straightforward  fashion 
as: 

N-3 

Step  1,  ( l (2)  (3)  (3))  + (2)  (2)  (2) 

i“l 

N-3 

Step  2,  ( l (2)  (2)  (3)  (N-i) ) + (2)  (2)  (2)  (2) 

i-1 

2 

which  totals  6N  +12N-66.  By  contrast,  for  the  decomposition  of  figure  2, 

= {i+l,i+2, .. . ,n)  which  is  exactly  the  same  as  if  the  network  was  fully 
connected,  and  it  follows  immediately  that  the  NXN  algorithm  requires  2N3 
operations.  This  example  makes  it  clear  that  the  choice  of  decomposition 
can  have  a profound  effect  of  the  efficiency  of  the  algorithm. 

For  an  arbitrary  network,  finding  the  optimal  decomposition  in  the 
sense  of  minimizing  the  required  number  of  operations  for  the  NXN  algorithm 
is  not  a trivial  problem  and  probably  can  only  be  solved  by  exhaustive 
comparision.  The  method  of  choosing  the  decomposition  for  the  examples 
which  are  presented  later  in  Section  IV  deviated  only  slightly  from  the 
following  heuristic  procedure: 
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1)  Label  a node  "1"  such  that  the  cardinality  of  C 
is  minimized. 

2)  For  every  ie{2,3,...,N>  in  order,  do  step  as 

a)  Given  the  nodes  which  have  been  labeled  "1" ,"2", . . . , 

"i-1",  label  an  unlabeled  node  "i"  such  that  the 
cardinality  of  is  minimized. 

The  effort  in  finding  the  decomposition  via  the  above  procedure  is  on 
the  same  order  as  doing  a shortest  route  computation  via  Floyd's  algo- 
rithm, and  as  a consequence  computer  time  savings  are  realized  only  when 
the  NXN  algorithm  is  iterated  several  times  for  the  same  topology. 

There  are  a large  number  of  networks  such  that  the  computation  time 
does  not  vary  widely  with  the  decomposition.  Such  networks  could  be  termed 
"locally  connected"  and  have  the  property  that  the  nodes  to  which  there  are 
direct  arcs  from  any  given  node  are  very  likely  to  have  direct  arcs  to  one 
another.  In  this  case,  the  nodes  could  be  numbered  very  rapidly  by  eye 
with  little  degradation  in  efficiency  (nodes  must  at  some  point  in  time 
be  assigned  a number  anyhow  in  order  to  communicate  the  topology  to  the 
computer) , and  in  the  first  shortest  route  computation  the  connection  sets 
could  be  established  with  very  little  effort.  In  f.ct,  the  only  modification 
to  the  NXN  algorithm  is  an  additional  step  which  is  included  just  before 
step  la: 

aa)  Initially  = 0;  for  je{i+l, i+2, . . . ,n} , do  step  bb: 
bb)  Include  j in  C . if  d.  . < °°. 

i I] 


The  additional  operations  required  by  this  step  number  ^N  which  is  quite 
modest  for  the  potential  gains. 
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Section  IV  The  IHU  Algorithm 


The  presentation  of  the  IHU  algorithm  requires  some  additional  defini- 
tions. For  this  algorithm,  a network  decomposition  is  defined  as  a division 
of  the  network's  nodes  into  ordered  subsets  S , S2,...,S^  such  that  for 

every  ieS  and  j£S0 , d. . = 00  if  |m-£|  > 1.  Every  node  of  the  network 
m 

belongs  to  exactly  one  subset.  The  submatrix  D contains  all  the  distances 

S . o 
1 m 

of  arcs  from  elements  of  S . to  elements  of  S and  has  dimension  Is  I x Is  I 

1 m ' i1  1 m1 

(where  |s|  means  the  cardinality  of  set  S) . Evidently,  D has  no  finite 

O . b 

1 m 

entries  in  the  case  |i-m|  > 1 (see  figure  3). 

Various  matrix  operations  will  be  performed  on  the  submatrices  to  gen- 
erate the  desired  shortest  distance  matrix.  Let  D *-  £ D mean  D 

0.0.  o.o,  0.0, 

11  11  11 

is  replaced  by  the  shortest  distance  matrix  computed  from  the  submatrix 

D . Define  A*B  = (min(a.  +b  .))  and  min(A,B)  = (min (a . . ,b . . ) ] . Also  let 
s.s.  1m  ml  11  ii 

11  m 

S,US_U...US  = (2  , and  if  m = k (where  k is  the  number  of  ordered  sets) 
12  mm 

then  S2  = £2,  = f2.  Define  the  conditional  shortest  distance  submatrix, 
m k 

D*  ((2  ) , as  the  shortest  distance  submatrix  under  the  restriction  that 
s.s.  m 
1 1 

all  the  intermediate  nodes  on  the  respective  conditional  shortest  routes 

are  members  of  f2  . In  the  case  n = k,  D*  _ ((2  ) = D*  ((2)  = D*  _ . 

m 0.0.  m s.s.  S.S. 

11  1 3 13 

Under  the  assumption  that  an  allowable  decomposition  has  been  given, 
the  following  algorithm  generates  all  the  shortest  distances  in  the  network 
(the  parenthetical  equality  to  the  right  of  each  step  is  the  claim  of  what 


I 


each  step  accomplishes) s 


Fiqure  3.  The  form  of  the  D matrix  for  the  IHU  algorithm  in 


the  case  k = 5.  If  the  decomposition  is  to  be 
acceptable,  the  shaded  submatrices  have  no  finite 
entries  prior  to  the  algorithmic  operations  on 
the  D matrix. 


! 


h 
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r.i 

k\ 


1)  D +■  E D 

S1S1  S1S1 


<-  °S  S «V  > 
11 


2)  For  every  ie{l,2, . . . ,k-l}in  order,  do  steps  a,  b,  c,  and  d: 


a) 

D +■ 

c c 

D *D 

c c c c 

(- 

Ds  s 

S . . b . 

1+1  1 

si+lsi  sisi 

Si+lSi 

b) 

D +• 

sisi+l 

D *D 

s.s.  s.s.A1 

11  1 1+1 

( = 

D* 

SiSi+l 

c) 

min  (D 

9 

D_  _ *1 

Si+lSi+l 

Si+lSi+l 

Si+lSi 

(B  )) 

Ds.s  > 

1 1+1 

(=  D*  {Q)) 

si+lsi+l  1 


d) 

D.  - ^ 5 D 

Ji+1  i+1  i+1  l+l 

(=  Ds  S ) 

Si+lSi+l  2+1 

3) 

For 

every  ie{k,k-l, . . . , 3, 2} 

in  order,  do  steps  a,  b,  and  c: 

a) 

D +■  D *D 

S.S.  , S.S.  S.S.  , 

1 1-1  11  1 1-1 

(=  D*  ) 

SiSi-l 

b) 

D + D *D 

S.  ,S.  S.S.  . S.S. 

l-l  l l l-l  l 1 

<=DS.  ,s.> 
1-1  1 

c) 

Ds.s.  , "min  (Ds.  lS.'D 
1 1-1  1-1  1 

s.  s Ds.s.  * (=Ds.  s } 

1-1  1 1 1-1  1-1  1-1 

4) 

For 

every  re{2, 3, . . . ,k-l}  in 

order,  do  step  a: 

a) 

For  every  i,je{l,2,...k} 

if  | i— j | = r,  do  step  b: 

b)  D„  „ +■  _ *D  _ 

s.s.  S.S  s S . 

13  1 P P 3 

(=Ds.s.> 
1 3 

where  p is  an  element  of  the  set  Q = {s+l,s+2, . . . , t-2,t-2}  for  s = min(i,j) 

and  t = max  (i,j)  such  that  |s  | < Is  I for  every  mGQ. 

p — m 

A rigorous  proof  of  the  algorithm  would  be  very  lenqthy  and  repetitious, 

and  the  interested  reader  is  referred  to  Hu's  work  [6]  for  exposition  of  a 

similar  proof.  Steps  1 and  2 are  bootstrapping  successive  diagonal  and 

first  off-diagonal  submatrices,  so  that  at  the  end  of  step  2,  D = D*  . 

SkSk  SkSk 


L . 
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Step  3 is  essentially  a backwards  form  of  step  2 and  replaces  the  diagonal 
and  first  off-diagonal  submatrices  with  the  respective  unconditional  shortest 
distance  submatrices.  Step  4 is  one  method  for  finding  the  unconditional 
shortest  distance  submatrices  corresponding  to  decomposition  sets  which  are 
separated  by  at  least  one  intermediate  set.  The  ordering  in  step  4 allows 
p to  be  any  element  of  the  set  Q,  and  the  particular  choice  of  p minimizes 
the  number  of  operations. 

If  one  assumes  that  the  shortest  distance  calculations  for  submatrices 
are  done  via  Floyd's  method  (requiring  2p3  operations  for  a p x p submatrix) 
and  that  the  pseudo-multiplications  are  done  in  a straightforward  manner 
(requiring  2pqr  operations  to  calculate  A»B  where  A is  dimension  p x q and 
B is  q x r) , then  the  number  of  operations  required  by  the  IHU  algorithm  is: 


Step  1, 
Step  2, 

Step  3, 
Step  4, 


2^1 1 3 

2 X l2|Si+li  ,Si 1 2 + |Si'  |Si*l|2  * |Si+l|3’ 

1k1 

2 l <2|S.|2  Is.^l  + |si_1|2|si|) 

2 E lSil  lS-jl  lSpl 

i»  j 

such  that  | i- j | > 1 


The  total  number  of  operations  is  then 


k-1 

2 ( l \S  US 
i=l 


1+1|3  - T Nl3  * I IsJ  Is  | Is  |) 

i=2  i,j 

such  that  | i- j | > 1 


- y,  ',i  riitdiWitottin  i ' Vir-f - ' 


One  may  compare  the  IHU  algorithm  to  other  versions  of  Hu's  algorithm. 

For  any  given  decomposition,  the  IHU  algorithm  requires  fewer  operations 

than  the  fastest  version  of  Hu's  algorithm  known  to  the  author,  which  is 

that  due  to  Yen  [11] . For  purposes  of  comparision,  an  example  which 

commonly  appears  in  the  literature  [5,6,11]  is  presented.  Let  |s^|  = 6 

for  i even  and  | | = t for  i odd.  Assume  6<_  t,  and  let  k,  the  number  of 

k+1 

sets,  be  odd.  Define  n = — — . In  this  case,  the  new  algorithm  requires 

3 2 2 2 2 2 3 

2(mt  + (m  +5m-6) t 6 + (2m  +2m-6) t6  + (m  -4m+5)6  ) operations.  Yen's 

3 2 2 2 2 2 3 

modification  requires  2(mt  + (m  +6m-7) t 6 + (2m  +10m-20)t6  + (m  +6m-14) 6 ) 

The  new  algorithm  is  faster  for  the  entire  range  of  interest,  i.e. 

t ^ 6 1 and  m ^ 2.  As  a particular  case,  let  5=  t and  m = 3|  the  IHU 

algorithm  requires  82t3  operations.  Yen’s  modification  requires  128t^ 

operations,  and  Floyd's  algorithm  requires  250t3  operations. 
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Section  VI  Some  Examples  Using  the  IHU  and  NXN  Algorithms 


In  this  section  several  examples  are  given  which  provide  insight  into 
the  classes  of  networks  for  which  the  NXN  and  IHU  algorithms  can  sub- 
stantially reduce  shortest  route  computation  time.  Although  no  examples 
are  presented  for  which  the  IHU  algorithm  is  faster  than  the  NXN  algo- 
rithm, they  do  exist.  Such  networks  form  a rather  small  and  special  class 
of  networks,  and  typically  may  be  decomposed  in  such  a manner  as  to  be  a 


variation  on  the  following  themes  |S^|  for  i odd  is  large  compared  to  |s^| 
for  i even,  and  if  jeS^  and  keS^  then  j and  k are  very  likely  to  have  direct 
arcs  to  one  another. 

The  first  example  is  an  old  version  of  the  ARPA  net  which  is  shown  in 
figure  4.  In  that  figure,  the  NXN  decomposition  is  defined  by  the  numbering 
of  the  nodes,  and  the  IHU  decomposition  is  defined  by  the  partitioning  of 
the  nodes  with  broken  lines.  This  network  lends  itself  to  NXN  decomposition 
due  to  the  high  number  of  nodes  which  have  arcs  directly  to  only  two  other 
nodes — a fact  which  keeps  the  cardinality  of  connection  sets  very  low. 

The  second  example  is  the  47  node  symmetric  network  shown  in  figure  5. 
This  network  is  not  "locally  connected"  to  a very  high  degree,  but  still 
the  NXN  algorithm  is  (perhaps  surprisingly)  efficient. 


i 

\\ 
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The  final  example  is  the  64  node  network  displayed  in  figure  6.  The 
density  of  arcs  is  perhaps  greater  here  than  in  the  other  examples,  but  a 
high  degree  of  local  connectivity  promotes  the  efficiency  of  the  NXN 
algorithm. 

Efficiency  is  measured  with  Yen's  implementation  of  Dijkstra's  algo- 
rithm as  the  standard.  Theoretical  efficiency  refers  to  the  relative 
savings  in  the  number  of  operations  required  to  perform  a shortest  route 
computation.  The  computation  times  for  the  IBM  370-168  to  execute  the 
Fortran  programs  of  various  algorithms  were  noted,  and  relative  savings 
are  referred  to  as  the  measured  efficiency.  The  comparisions  of  the 
various  algorithms  in  performing  shortest  route  calculations  on  the  three 
sample  networks  are  summarized  in  tab.e  1.  The  Fortran  programs  were 
complied  by  the  IBM  G1  compiler;  and  each  algorithm  not  only  computed  the 
shortest  distance  matrix,  but  also  computed  a routing  matrix  which  specified 
the  next  node  from  each  node  on  the  shortest  route  to  any  other  node. 


Fiqure  5.  A 47  node  symmetric  network  (nodes  are  connected 

to  first  and  seventh  nearest  neiqhbors  by  arcs) . HXN 
decomposition  is  indicated  by  node  labeling.  IHU  decom- 
position sets:  S = {l},  S = {21,47,20,39},  S = (34, 

40,30,26,2,4,5,8)7  S = { 14 , 17 , 18 , 29, 22 , 23, 42 , 41 , 32 , 25 , 
36,43},  S = {46,33,15,37,31,38,44,3,6,7,9,29,11,12}, 
and  = 113,27,16,19,45,35,10,24}. 


1 

Number  of 
operations 

27040 

154630 

393216 

Dijkstra's 

shortest 

Theoretical 

efficiency 

1.00 

1.00 

1.00 

route 

algorithm 

Computation 
time  (seconds) 

.090 

.450 

1.115 

Measured 

efficiency 

1.00 

1.00 

1.00 

Number  of 
operations 

6948 

83880 

124722 

IH 

algorithm 

Theoretical 

efficiency 

3.89 

1.84 

3.15 

Computation 
time  (seconds) 

.020 

.195 

.290 

Measured 
ef ficiency 

4.50 

2.31 

3.84 

Number  of 
operations 

2828 

28608 

42416 

NX’: 

alqorithm 

Theoretical 

efficiency 

9.56 

5.41 

9.27 

Computation 
time  (seconds) 

.015 

.115 

.165 

y. 


Table  1 


Comparative  performance  of  three  different  shortest 
route  alqorithms  on  the  three  sample  networks. 
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Appendix 

Section  IA  Introduction 


This  appendix  describes  and  lists  the  program  which  provided  the 
computational  experience  cited  in  this  paper.  The  program  of  section  VI IA 
reads  the  topology  of  the  network,  finds  a decomposition  for  the  IHU  and 
NXN  algorithms,  solves  a sample  shortest  route  problem  via  each  algorithm 
and  the  Dijkstra  algorithm  in  order  to  compare  computation  times,  and 
calculates  the  number  of  operations  required  by  each.  Typically,  an 
application  of  these  programs  requires  at  most  two  of  the  listed  subroutines — 
one  to  decompose  the  network  and  one  to  calculate  all  the  shortest  routes. 

The  decomposition  subroutine  needs  to  be  called  only  one  time  for  any  given 
topology  since  a new  set  of  data  cards  are  punched  by  the  decomposition  sub- 
routines which  record  the  appropriate  decomposition  information.  In  this 
appendix,  a hybrid  notation  will  be  employed  which  is  a combination  of 
that  used  in  the  body  of  this  report  and  that  used  in  the  Fortran  programs. 

The  definitions  of  all  Fortran  terms  are  given  in  the  comment  cards  at  the 
beginning  of  the  program  listing  that  is  found  in  section  VIIA. 


I 
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Section  IIA  Bookkeeping  for  Shortest  Routes 

The  algorithms  which  are  listed  not  only  find  the  shortest  distances 
between  every  pair  of  nodes  in  the  network,  but  they  also  record  the 
shortest  routes.  The  method  which  is  used  for  this  purpose  is  establishing 
a "next  node"  matrix  where  HX(I,J)  is  the  next  node  on  the  shortest  path 
from  node  X to  node  J.  Initially,  NX(I,J)  = J for  every  existing  arc 


(I,J)  , and  every  time  the  operation,  d^  +■  min  (d^,  dijt+dxj^  is  Performed 
such  that  diJc+dkj  the  distinct  minimum,  then  the  algorithm  makes  the 
replacement  HX(I,J)  «-  NX(I,K) . For  the  remainder  of  this  appendix,  the 
algorithms  are  discussed  only  in  terms  of  the  shortest  distance  problem. 


Section  IIIA  The  Main  Program 


The  main  program  reads  in  the  topology,  assigns  arc  numbers  and 


provides  the  control  for  its  specific  purpose,  i.e.  to  compare  the 


2 lists  the  data  cards  which  communicate  the  topology  of  the  network  to 


the  program.  The  first  card  is  a header  which  provides  the  name  of  the 


network  and  the  values  for  NH,  MXHU,  MNXIJ,  MAXPRI  and  NFORBD.  The  second 


card  says  that  node  ”1"  has  "2"  outgoing  arcs  which  terminate  on  nodes 


There  is  one  such  card  for  each  node  in  succession 


7 NODE,  13  ARC  EXAMPLE  NT  7 2 2 7 0 


Table  2.  Topology  Cards  For  The  Network 

of  Figure  7 
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Section  IVA  Subroutine  DIJKST 


Subroutine  DIJKST  is  an  implementation  of  Dijkstra's  algorithm  sug- 


gested by  Yen  [12].  The  algorithm  can  be  floated  to  perform  the  operation. 


D *■  £D  , in  the  case  that  if  j is  a node  number  such  that 

s , s , s . s . 

i i 11 


min  k <_  j <_ max  k,  then  jes..  When  the  call  to  the  subroutine  DIJKST  is 
kes . kes . 1 

made,  for  this  case,  then  NB  = min  k and  NF  = max  k.  If  the  operation, 

kes . kes 

D D*,  is  to  be  performed  via  Dijicstra's  algorithm,  then  NB  = 1 and 


NF  = NN. 
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Section  VA  Subroutines  DECIHU  and  IHU 

Subroutine  DECIHU  decomposes  the  network  for  the  IHU  algorithm  which 
is  implementated  in  subroutine  IHU.  The  method  of  decomposition  is  that 
of  Section  V.  Figure  8 shows  the  network  as  decomposed  by  DECIHU  with 
the  new  node  numbers  as  printed  out.  Table  3 shows  the  cards  punched  by 
DECIHU  which  record  the  decomposition  information  and  describe  the  topology 
in  terms  of  the  new  node  numbers.  Again,  the  first  card  is  a header  with 
the  title  of  the  network,  a "1"  which  says  the  cards  were  punched  by  DECIHU 
and  a "3"  which  is  the  number  of  IHU  sets.  The  second  card  says  that  node 
"1"  has  "2"  outgoing  arcs,  is  a member  of  set  number  "1"  (the  next  two 
zeros  have  no  significance) , and  the  outgoing  nodes  are  to  nodes  "2"  and 
‘’3”;  and  so  forth.  The  ninth  card  is  a header  for  NTWIXT  which  starts 
on  the  next  card.  From  them,  NTWIXT(1,1)  = "0",  NTWIXT (1,2)  = "0", 

NTWIXT (1, 3)  = "2",  NTWIXT (2,1)  = "0",  etc.  The  information  on  these  cards 
define  the  variables  found  in  the  common  block  IHUSTF,  and  these  values 
are  given  in  Table  4. 

Subroutine  IHU  is  a straightforward  implementation  of  the  IHU  algorithm 

as  presented  in  Section  IV.  The  operations,  Dg  *■  £Dg  , are  performed 

i i i i 

via  subroutine  DIJKST. 


i 


■ - - - 
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Figure  8.  Node  Renumbering  and  Partitioning  by 
Subroutine  DECIHU  for  the  Network  of 
Fiqure  7. 
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7 NODE,  13  ARC  EXAMPLE  NT  7 1 3 


1 2 1 0 0 2 3 


2  2 2 0 0 1 6 


3 2 2 0 0 6 5 


4 2 2 0 0 1 7 


5 2 3 0 0 4 7 


6 1 3 0 0 5 


7 2 3 0 0 2 4 


NTWIXT  FOR  7 NODE, 13  ARC  EXAMPLE  NT 


00  2 0 0030  0 


Table  3.  Cards  punched  by  subroutine  DECIHU  which  relate 

the  IHU  decomposition  information  and  the  topology 
in  terms  of  the  new  node  numbers  for  the  network 
of  Figure  7. 


Nl(l)  = 1 


Nl(2)  = 2 


N1 ( 3)  = 5 


N2(l)  = 1 


NTWIXT (1,1) 
NTWIXT (2 , 1) 


N 2(2)  » 4 
NTWIXT ( 2,2) 
NTWIXT (2, 2) 


N2 (3)  = 7 
NTWIXT (1,3) 


NTWIXT  (2,3) 


NTWIXT  0,1) 


NTWIXT(3,2) 


NTWIXT(3, 3) 


NS  = 3 


Table  4.  Values  of  the  variables  in  labeled  common 

block  IHUSTF  which  may  be  deduced  from  cards 
in  Table  3 for  the  network  of  Figure  8. 


Subroutine  NXN  is  a general  implementation  of  the  NXH  algorithm  for 


i 


the  case  in  which  all  the  arcs  in  the  network  are  not  necesarily  duplex. 
Two  connection  sets  are  established  for  each  node — one  for  outgoing  con- 
nections and  ore  for  incoming  connections.  Define  C°  as  the  outgoing 
connection  set,  i.e.  jeC°  if  there  exists  a path  P from  i to  j such  that 
C (P)  < 00  and  every  intermediate  node  k satisfies  k < i.  Similarly, 

define  C1  as  the  ith  incoming  connection  set.  The  NXN  algorithm  takes 
1 

this  form: 

1)  For  every  ic{l, 2, . . . ,NN-2}  in  order,  do  step  a: 

I o 

a)  For  every  je  and  ke  C^,  do  step  b: 

b)  d.,  min  (d  , d..+d  ) 

3k  3k  31  lk 

2)  For  every  ie{NN-2,NN-3, . . . ,l}  in  order,  do  step  a: 

a)  For  every  je{i+l,  i+2,...,NN>  , do  steps  b and  c: 

b)  For  every  meC°,  d,  . min  (d.  +d  d.  .) 

l 13  1m  m3  13 

c)  For  every  keC1,  d . . min  (d..,  d.,+d,  .) 

1 31  31  3k  ki 

The  method  DECNXN  uses  for  decomposing  the  network  is  given  in  Section 

III  with  the  alteration  that  nodes  are  chosen  in  order  to  sucessively 

minimize  |c?|  + |c*|.  For  the  network  of  Figure  7,  the  new  node  number- 
1 ' 1 

ing  which  implies  the  decomposition  is  shown  in  Figure  9.  The  cards 
punched  by  DECNXN  which  contain  topology  information  in  terms  of  new  node 
numbers  and  the  decomposition  information  are  shown  in  Table  5. 


j 


I 
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The  interpretation  of  the  cards  is  now  more  difficult  but  should  be 
clear  by  the  program  in  Table  6 which  reads  in  the  cards  of  Table  5,  sets 
up  arc  numbers,  and  prepares  the  decomposition  information  for  DECNXN. 

One  feature  of  the  program  not  yet  discussed  is  that  of  NFORBD  which 
is  an  input  variable.  If  a network  is  "locally  connected"  except  for  a 
few  nodes,  they  should  be  numbered  last  and  suppressed  from  being  assigned 
new  node  numbers  which  are  low  by  establishing  NFORBD  as  the  cardinality 
of  the  set  of  such  nodes. 


7 NODE,  13  ARC  EXAMPLE  NT  7 2 

12  32  323  2 5 3 

21  21  355  4 3 

32  22  341  5 4 6 

42  21  332  56  7 

52  12  267  6 7 

6 2 1 0 0 3 7 

7 2 1 0 0 4 6 

Table  5.  Cards  punched  by  subroutine  DECNXN  for 

the  network  of  Figure  7 which  contain  decom- 
position information  for  the  NXN  algorithm 
and  topology  information  in  terms  of  the  new 
node  numbers. 
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Section  VIIA  Program  Listing 

The  program  of  this  section  provided  the  computational  results  of 
this  paper.  The  program  is  generously  commented  and  should  be  transparent 
when  studied  along  with  this  appendix.  In  general,  clarity  was  sacrificed 


for  speed  only  in  the  subroutines  DIJKST,  NXN,  and  IHU. 
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DO  3304  K=INDEX1 , LV 
LC=NC  (K) 
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